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ABSTRACT 

We present a high angular resolution map of 850 /im continuum emission of 
the Orion Molecular Cloud-3 (OMC 3) obtained with the Submillimeter Arraj0 
(SMA); the map is a mosaic of 85 pointings covering an approximate area of 
6'.5x2'.0 (0.88x0.27 pc). We detect 12 spatially resolved continuum sources, 
each with an H2 mass between 0.3-5.7 M and a projected source size between 
1400-8200 AU. All the detected sources are on the filamentary main ridge 
(tih 2 >10 6 cm -3 ), and analysis based on the Jeans theorem suggests that they 
are most likely gravitationally unstable. Comparison of multi-wavelength data 
sets indicates that of the continuum sources, 6/12 (50 %) are associated with 
molecular outflows, 8/12 (67 %) are associated with infrared sources, and 3/12 
(25 %) are associated with ionized jets. The evolutionary status of these sources 
ranges from prestellar cores to protostar phase, confirming that OMC-3 is an 
active region with ongoing embedded star-formation. We detect quasi-periodical 
separations between the OMC-3 sources of ~17"/0.035 pc. This spatial distri- 
bution is part of a large hierarchical structure, that also includes fragmentation 
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scales of GMC (^35 pc), large-scale clumps (?^1.3 pc), and small-scale clumps 
(~0.3 pc), suggesting that hierarchical fragmentation operates within the Orion 
A molecular cloud. The fragmentation spacings are roughly consistent with the 
thermal fragmentation length in large-scale clumps, while for small-scale cores it 
is smaller than the local fragmentation length. These smaller spacings observed 
with the SMA can be explained by either a helical magnetic field, cloud rota- 
tion, or/and global filament collapse. Finally, possible evidence for sequential 
fragmentation is suggested in the northern part of the OMC-3 filament. 

Subject headings: ISM: clouds — ISM: individual (OMC-3) — stars: formation 
- Stars: protostars — ISM: structure — ISM: individual (Orion A) 

1. INTRODUCTION 

The key steps of the star formation process such as large-scale gas infall, disk accretion, 
and mass ejection have been studied mainly in nearby isolated star-forming regions by vari- 
ous methods such as the analysis of: Spectral Energy Distributions (SEDs), gas kinematics, 
density structures of parental cores, and molecular outflows/jets (e.g., Evans et al. 1999; 
Bachiller & Tafalla 1999; Myers et. al. 2000 and therein). However, infrared studies suggest 
that the majority of stars (80-90%) form as members of embedded clusters (N*>35) within 
nearby giant molecular clouds (GMCs; located at d<2 kpc), rather than in isolated star 
formation environments (e.g., Lada & Lada 2003; Porras et al. 2003). Hence, studies of 
the formation and evolutionary processes of embedded clusters, as well as associated GMCs, 
are crucial to understand the general features of star formation. Moreover, these embedded 
clusters, identified as a group of infrared excess sources, are often associated with dense 
filamentary structures (i.e., N# 2 >10 22 cm" 2 ), and most likely represent the youngest asso- 
ciation of the embedded stellar-clusters (Myers 2009). In order to observationally investigate 
the initial conditions of embedded clusters, it is essential to spatially resolve the individual 
star-forming cores within molecular filaments, obtain kinematic information, and investigate 
the fragmentation properties and physical conditions of each of the cluster members. 

The Orion A molecular cloud (d = 414 pc; Sandstrom et al. 2007; Menten et al. 2007; 
Hirota et al. 2007; Kim et al. 2008), which contains a large amount of gas (~5xl0 4 M ), 
is the nearest GMC, and one of the best-studied embedded cluster-forming regions at all 
observable wavelengths (e.g., Bally et al. 1987; Allen et al 2007; Tatematsu et al. 1993, 
Johnstone & Bally. 1999; Tsujimoto et al. 2005; O'Dell et al. 2008). A large number 
of molecular line observations have been performed: in 12 CO by Maddalena et al. (1986), 
Castets et al. (1990), in 13 CO by Bally et al. (1987); Buckle (2012), in C 18 by Dutrey 
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et al. (1991; 1993); Buckle (2012), in CS by Tatematsu et al. (1993), in NH 3 by Castets 
et al. (1990), Harju et al. (1991), Cesaroni & Wilson (1994), and Wiseman & Ho (1998), 
and in H 13 CO+ by Aso et al. (2000), Ikeda et al. (2007), Tatematsu et al. (2010). These 
observations revealed the distribution of the dense gas within the Orion A GMC down to 0.1 
pc scale. In addition to molecular line observations, (sub) millimeter continuum observations 
performed by single-dish telescopes have revealed a remarkable filamentary structure over 
the entire Orion A molecular cloud with a length of ~7 pc. A large number of dust contin- 
uum sources, candidate Class 0/1 sources, were detected, and properties of the dust emission 
along the Orion Molecular Cloud (OMC) filaments were investigated (Chini et al. 1997; Lis 
et al. 1998; Johnstone & Bally 1999; Nielbock et al. 2003; Nutter & Ward-Thompson 2007). 
Furthermore, near- to mid- infrared surveys by Spitzer MIPS/IRAC showed that a num- 
ber of infrared excess sources (YSO candidates) are actually located on the OMC filaments 
with a wide stellar mass range, from brown dwarfs to massive stars (Nielbock et al. 2003; 
Allen et al. 2007; Peterson et al. 2008; Peterson k Megeath 2008; Sadavoy et al. 2010; 
Megeath et al. 2012). Moreover, recent high-angular resolution radio observations, targeted 
at relatively bright (sub) millimeter sources (i.e., intermediate-mass protostellar candidates), 
directly resolved detailed spatial structure of protostellar cores and associated molecular 
outflows in the OMC filaments. The detected sources showed a variety of properties such 
as different structures, kinematics, and evolutionary phases (e.g., Williams et al. 2003; Za- 
pata et al. 2004; 2005; 2006; 2007; Matthews et al. 2005; Takahashi 2006; 2008a; 2008b; 
2009; 2012; Shimajiri et al. 2008; 2009; Takahashi & Ho 2012). However, so far only a 
very limited number of millimeter and submillimeter observational studies have been made 
which achieve an angular resolution comparable to, or better than, existing IR data (i.e., 
<a few arcsec resolution, which can spatially resolve <1000 AU scale structures at the dis- 
tance of Orion). Hence, the physical relation between IR sources and the small-scale gas 
structure (D< a few 1000 AU and high-density gas: n H2 >W 5 cm -3 ), which can be spatially 
resolved with the millimeter /submillimeter interferometric observations, remains unclear. It 
is also important to spatially resolve individual star forming sites, or cores, in order to study 
the (IMF precursor) core mass function (CMF). Analysis of the CMF could lead to the un- 
derstanding of the similarities and differences in the IMF of entire stellar-clusters which may 
have been imposed at the initial evolutionary stage. 

In order to observationally study the nature of embedded sources and their associated 
filaments, as a first step, we report 850 fim continuum observations of the molecular filaments 
toward the OMC-3 region located in the northern part of Orion A GMC with the largest 
mosaic observations ever published using the SMA (6'.5x2'.0; 0.88x0.27 pc). The SMA, 
equipped with 345 GHz receivers is one of the most powerful facilities to spatially resolve 
prestellar /protostellar cores and circumstellar disks through the thermal dust continuum 
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emission. The goals of the paper are to (i) spatially resolve the individual star-forming 
cores with the angular resolution of 4". 5 (~1900 AU) and identify the prestellar and proto- 
stellar candidates with a mass sensitivity of 0.15 M , (ii) compare the physical properties 
of (sub)millimeter sources with existing multi-wavelength data and study the evolutionary 
status of each identified source, and (iii) quantify the physical properties of the filaments 
such as fragmentation scale lengths by comparing our SMA images with previous single-dish 
images. 

We explain in Section 2 the strategies of the observations, data reduction, and imaging. 
We present in Section 3 the synthesized images and the derived physical parameters of the 
filaments. In Section 4, we discuss the physical properties and the evolutionary status of 
identified cores and the derived fragmentation properties of the OMC filaments. Finally, we 
summarize the paper and discuss future prospects in Section 5. 

This paper serves as a starting point for a series of high-angular resolution OMC stud- 
ies focused on the filamentary structures and their embedded sources. Our studies, which 
will cover the entire system of OMC filaments in the continuum, as well as molecular line 
emissions in the (sub) millimeter (described in Section 2.1), will include a statistically signifi- 
cant number of embedded sources. Once the Atacama Large Millimeter /submillimeter Array 
(ALMA) and the Jansky Very Large Array (JVLA) are fully operational, the sensitivity and 
resolution of their data will be greatly improved so that the wide-field mosaic observations 
will be a common approach, especially for GMCs and nearby extragalactic objects. Hence 
our studies would be very important in developing a strategy to investigate structures of 
molecular clouds and associated proto-clusters on pc-scale down to AU scale. 



2. OBSERVATIONS AND DATA REDUCTION 

2.1. SMA Observations 

The observations were carried out at 351 GHz using the eight 6 m antennas of the SMA 
(Ho et al. 2004) in the compact configuration on January 2 and October 21, 2008 and in the 
subcompact configuration on August 29, 2008. The primary beam at this frequency is 36", 
and the 85 pointing centers were arranged with sampling in a nyquist grid spaced by 18" with 
a position angle of 45° (Figure 1). We observed each field for one minute during each of 4 - 18 
cycles through all the fields, giving a total on-source integration time of 4-18 minutes. Both 
the LSB (z/ c =346 GHz) and USB (z/ c =356 GHz) data were obtained simultaneously with the 
90° phase switching technique by the digital spectral correlator, which had a bandwidth of 
2 GHz for each sideband. 
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The SMA correlator covers 2 GHz bandwidth in each of the two side-bands separated 
by 10 GHz. Each band is divided into 24 chunks of 104 MHz width, which were covered 
by fine spectral resolution (256 channels correspond to a velocity resolution of 0.36 km s -1 
in our observing setting). In addition to continuum observations, molecular lines such as 
HCO+(4-3), H 13 CO+(4-3), CO (3-2), SiO(8-7), and several CH 3 OH and SO transitions were 
simultaneously obtained with the 2 GHz band width in each sideband. These molecular line 
data will be presented in a future article. 

The continuum data from both sidebands were combined to obtain a higher signal-to-noise 
ratio. The combined configurations of the arrays provided projected baselines ranging from 
9.7 to 88 kA. Our observations were insensitive to structures larger than 17" (7000 AU) at 
the 10% level (Wilner & Welch 1994). Passband calibration was achieved by observing the 
quasars, 3C 273 and 3C 279. Amplitudes and phases were calibrated by observations of 
0530+135 and 0423-013, whose flux densities were determined relative to Titan and Nep- 
tune. The overall flux uncertainty was estimated to be ~20%. The pointing accuracy of the 
SMA observations was ~3". The raw data were calibrated using MIR, originally developed 
for the Owens Valley Radio Observatory (Scoville et al. 1993) and adopted for the SMA. 
The SMA observational log for each date and the observational parameters are summarized 
in Table 1 and Table 2. 



2.2. Mosaic Imaging Procedures 

In order to produce the mosaic continuum image, certain frequencies (which include 
possible molecular lines described above) are subtracted from the visibility data. The effec- 
tive bandwidth for the continuum emission is approximately 3.2 GHz. CLEANed images in 
each pointing were made using DIFMAP with a Gaussian taper (HWHM=25 kA). Due to 
residual phase errors, it was difficult to remove the increased noise, which is greater than 5a 
level, around the four brightest compact sources, SMM 2, SMM 3, SMM 8, and SMM 11. 
These enhanced noise structures are due to the limited dynamic range which can be achieved 
against the strong peak signal levels of those bright sources. In order to improve the dynamic 
range of the images, we reduced the phase noise errors, by applying a self-calibration pro- 
cedure toward the several fields which include either SMM 2, SMM 3, SMM 8, or SMM 11 
within their primary beams. Within the integration time of 4-18 min., the theoretical noise 
levels were estimated to be 13- 25 mJy beam" 1 . After self-calibration, the theoretical noise 
levels were achieved around SMM 2, SMM 3, and SMM 11. Images around SMM 8 (i.e., 
nine mosaic points centered on SMM 8) achieved 1.1 to 2.0 times theoretical noise levels. 
This corresponds to a dynamic range of ~200 around SMM 8. 



After the CLEANing and self-calibration processes, the 85 images were imported into 
Miriad (Wright & Sault 1993) for primary beam corrections, and then combined into a 
single image using the miriad task "linmos" . Since each image had slightly different angular 
resolution, all the images were convolved to the final beam size of 4.5" FWHM (corresponding 
to 1900 AU). The achieved rms noise level for the entire mosaic image was estimated to be 
20 mJy on average. 



3. RESULTS 

3.1. Identification of SMA 850 |im Compact Continuum Sources 

Figure 2 shows the SMA 850 /zm continuum in the left panel, and in the right panel we 
compare the SMA image with the previously published JCMT/SCUBE 850 /im image by 
Johnstone & Bally (1999). We clearly detected 12 compact SMA continuum sources embed- 
ded within the filament with a signal-to-noise ratio of greater than 5 (flux more than 100 mJy beam 
denoted by white contours in the right panel of Figure 2). All the detected continuum sources 
are found along the main ridge of the OMC-3 filament, with no detections off the main ridge. 
This indicates that the relatively massive sources in this region are likely formed only within 
the densest part of the filaments. In order to measure the position, peak flux, total flux 
density, and deconvolved size of the continuum source, 2D Gaussian fitting in AIPS (task 
"JMFIT") was used. The measured results and their uncertainties are listed in Table 3. The 
spatial structures of the individual sources are presented in Figure 3. 

Eight out of the twelve SMA continuum peaks (SMM 1-6, SMM 7, and SMM 10) coin- 
cide with (sub) millimeter continuum sources previously detected with the JCMT/SCUBA, 
CSO/SHARC, and IRAM/MPIfR, to within the single-dish positional accuracies of ~±5" 
(Chini et al. 1997, Lis et al. 1998, Johnstone & Bally 1999). SMM 10 was detected in the 
previous 3.3 mm interferometric observations at the 3a level as MMS 7-NE (Figure 2a of 
Takahashi et al. 2006). Figure 1 of Takahashi et al. (2009) shows an image where SMM 
7 (MMS 6-NE) is spatially resolved using SMA compact configuration data with robust 
weighting. The sources, SMM 9 and SMM 12, were detected for the first time with our SMA 
850 fim continuum observations. For the other eight sources, this is the first time that they 
have been spatially resolved and positioned accurately. On the assumption that the 850 /im 
emission comes from optically thin dust emission and that the temperature distribution of 
the dust continuum emission is uniform, we can estimate the source mass by 
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where k\ is the mass- absorption coefficient of dust grains, B\(T) is the Planck function 
for a dust temperature T^ ust , Fx is the total flux density of the continuum emission, and 
d is the distance to the source (414 pc; Menten et al. 2007). Adopting a dust opac- 
ity of /t A =0.037 cm 2 g- 1 (A/400 /im) ~? (Keene et al. 1982), (5 = 1.6 (Johnstone & Bally 
1999), a dust temperature of Td us t = 20 K (Cesaroni & Wilson 1994, Chini et al. 1997; 
Tatematsu et al. 2008), a gas-to-dust ratio=100, and measured total flux densities (Table 
3), the masses of the continuum sources are estimated to range from 0.3 to 5.7 M Q . High 
mean molecular hydrogen volume densities were also derived toward these sources with a 
range between 3.2xl0 6 < nu 2 <3.3xl0 8 cm~ 3 under the assumption of a spherical geometry 
(nn 2 = Mu 2 / (^irr 3 [xmii)) . Here, /i and m# are the hydrogen mean molecular weight (2.33) 
and the hydrogen mass, respectively, and r is the radius of the source, taken as the geometric 
mean of the major and minor axes of its deconvolved size. The derived physical parameters 
of these sources are summarized in Table 4. 

Note that the total flux of the 850 /im continuum emission might include contamination 
from ionized emission, which comes from free-free continuum jets of protostars. From 3.6 cm 
observations (Reipurth et al. 1999), and with the assumption that the typical spectral index 
of free-free jets is 0.6 (Anglada et al. 1998; Reynolds 1986), we find that this contamination 
at 850 /im is negligible (i.e., <0.4 %) for all the detected OMC-3 sources. The continuum 
emission at 850 /im must therefore come from the thermal dust. 

In order to estimate the flux concentration factors, C(sma/jcmt), the SMA data were 
first convolved to the JCMT 14" beam, and then the peak intensity within the JCMT 14" 
beam (Jy beam -1 ) are used for deriving the ratios. The flux concentration factor for SMM 8 
(or MMS 6-main; Takahashi et al. 2009;2012b) is 86%, approximately twice higher than 
those of the other detected continuum sources (The mean and median value of the flux 
ratios were estimated to be 41% and 35 %, respectively). SMM 8 shows small deconvolved 
sizes compared with most of the other SMA sources and has the highest flux among the 12 
detected sources (Table 4). No clear difference of the flux concentration ratio was detected 
between the prestellar cores (i.e., core not associated with any outflows nor any IR sources; 
Table 5) and the protostellar cores (i.e., core associated with either outflows or IR sources; 
Table 5) except for the higher ratio value derived in SMM 8. The evolutionary stage of 
detected sources will be discussed in Section 4.2. 

Moreover, to determine the flux- recovering ratio for the entire filament, we proceeded 
in the following manner. First, we calculated the total SMA flux over the entire image (after 
convolving the image with the JCMT beam), F SM a (~16 Jy/14" beam); this corresponds to 
the sum of the fluxes of each of the SMA sources since we detected no extended emission 
from the filament itself in our SMA image. Second, we determined the JCMT contour level 
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that corresponds to the the SMA 5a detection level (100mJy/4".5 beam), i.e., 0.97Jy/14" 
beam; this was used as a lower threshold for extended emission in the JCMT image. We 
then summed up the flux within this contour in the JCMT image, -Fjcmt (~103 Jy/14" 
beam). Finally, we divided the two fluxes, -Fsma/-Pjcmt, and found that only 16% of the 
total JCMT flux was recovered in our SMA data. 

It is technically possible to combine the SMA 850 /im data with the JCMT/SCUBA 
data in order to recover the flux from short antenna spacings. However, we did not apply 
this method due to the following reasons: (i) JCMT/SCUBA is a bolometer and covers 
a much wider frequency range as compared to the SMA (A5jcmt^7AI?sma), and (ii) 
JCMT/SCUBA frequency coverage includes several strong molecular lines such as CO, 
HCO + , and SO. Hence, the JCMT continuum emission map may be contaminated by emis- 
sion from these molecular lines. 



3.2. Physical Properties of the Continuum Sources 

In order to further discuss the properties of the detected sources, it is important to know 
the limits of our SMA observations. 

Figure 4(A) shows the mass distribution of the detected sources in OMC-3 (filled his- 
togram). The median and mean value of the source masses were estimated to be 1.6 M and 
1.9 M , respectively. Among the detected sources, SMM 8 (or MMS 6-main in Takahashi 
et al. 2009) is a source which has a gas mass 3.5 times greater than the median value. Note 
that SMM 8 has the largest mass, the largest flux concentration factor, and one of the most 
compact source sizes, indicating a very different nature as compared to the other sources in 
OMC-3. Detailed properties of this source are discussed in Takahashi et al. (2009); (2012), 
and Takahashi k Ho (2012). 

Since the number of detected sources were limited, and the error bars for the distribution 
scale as \fN (a Poisson distribution was adopted for the errors), it is difficult to discuss the 
characteristic mass from our present results. A full analysis of core census will be our future 
goal with a more extensive sample. However, an interesting question is whether we can 
extrapolate the mass distribution below the detection limit of our observations. According 
to our previous calculation, the JCMT "residual flux" is -Fjcmt — -^sma=87 Jy. This total 
flux corresponds to an H 2 mass of 98 M Q , with an assumption of /3=1.6 (Johnstone & Bally 
1999). This suggests that the residual gas has the potential to form numerous cores (~98 
of 0.1 M cores) even with a 10 % star formation efficiency (estimated from the nearby 
embedded clusters; e.g., Lada & Lada 2003). 
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We consider that it is possible to detect sources less than 0.1 M because recent SMA 
1.3 mm continuum observations have actually detected 31 circumstellar disks with a mass 
range from 0.0028 - 0.0662 M toward the same Orion star-forming region (28 in the Trapez- 
ium cluster and 3 outside of the Trapezium cluster; Mann & Williams 2009 a, b) as shown 
in Figure 4A (open histogram). The relatively shallow and large mosaic observations we 
performed, as compared to the targeted observations performed by Mann & Williams (2009 
a, b), make it difficult to detect these less massive circumstellar matter around the relatively 
evolved young stars (Class II sources). Moreover, note that these disks have much smaller 
sizes (<300 AU) than that of the sources detected in our SMA observations. 

Figure 4(B) shows the source size distribution. The source sizes are mainly distributed 
from 1000 to 4000 AU. In our OMC-3 sample, the minimum detectable source size is roughly 
limited by the present SMA angular resolution B However, as noted before, the results of the 
Orion disk survey suggest the presence of more compact sources (D < 300 AU; Vicente & 
Alves 2005; Mann & Williams 2009 a, b). The maximum detectable source size is determined 
by the shortest baseline in the database, which is A/-D m i n .~19" or 7800 AU. Further discus- 
sion on the detection of extended sources can be found in Wilner & Welch (1994). They 
used Gaussian models with limited uv sampling. If we assume Gaussian shaped sources, our 
observations are sensitive to ~17" (7000 AU) structures at the 10 % level. The source size 
of SMM 5, shown as the largest circular source size in Figure 4(B), is larger than the SMA 
detectable size - this is because we used a two-dimensional Gaussian fitting method for de- 
termining the source size. SMM 5's weaker peak intensity (compared to other sources), and 
its non-Gaussian structure, resulted in a two-dimensional Gaussian fit with a large FWHM. 

In the residual image, the mean flux is 1.4 Jy within the JCMT 14" beam. With the 
assumption of uniform flux distribution, the mean flux within the SMA 4". 5 beam is esti- 
mated to be 144 mJy. This value is larger than the SMA 5cr level (100 mJy), which is the 
threshold for source identification. However, unless this missing flux arises from residuals 
that are smaller than the maximum detectable size, it will not be detected. Therefore, the 
missing flux most likely arises from the more extended filamentary structure (>19"). 



2 The minimum detectable size is also related to the signal-to-noise ratio of the detected sources. Since 
enough signal-to-noise ratio (>50) has been achieved in our observations, we have succeeded to spatially 
resolve the structures less than the synthesized beam. 
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4. DISCUSSION 

4.1. Star formation in OMC-3 

As shown in the previous section, potential star forming cores have been spatially re- 
solved at 850 fxm with the SMA. The origin of the 850 /im continuum emission comes mostly 
from thermal dust emission, suggesting dense gas associated with circumstellar envelopes. 
Physical properties of the detected continuum sources and the star formation activities in 
the OMC-3 region are discussed in the following subsections. 



4-1.1. Gravitational Stability of the Continuum Sources 

In order to discuss the physical properties of the detected SMA continuum sources, their 
measured source sizes are plotted as a function of the estimated hydrogen number density. 
Jeans analysis can provide a useful way to discuss the gravitational stability of the observed 
sources. With the assumption of an infinite and homogenous medium, the Jeans length is 
described as follows (Jeans 1902), 



*frag. 



7TC? 



G Po 



(2) 



where G is the gravitational constant, p is the mean density, and c s is the sound speed 
(the relation between the sound speed and the temperature of the gas are described as 
where T gas is the gas temperature and k is the Boltzman constant. The 



Jeans length corresponds to the critical scale when the gas cloud becomes gravitationally 
unstable, i.e., if the source size is larger than the Jeans length (-D CO re>Af rag .), the core becomes 
gravitationally unstable and starts to collapse. The Jeans critical number density is described 
as a function of the Jeans length and gas temperature as: 

n H2 (A frag ,T gas ) = ^-^^- (3) 

In Figure 4(C), the Jeans critical number density is plotted as a function of source size 
for different gas temperatures (T=10 K, 20 K, and 30 K) and is overlayed with the observed 
results. Here, we consider that the gas temperature within OMC filaments is approximately 
20 K (e.g., Cesaroni & Wilson 1994). The plotted results clearly show that most SMA 
sources are located very close to the critical line (T gas =20 K) and several sources have a 
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larger size than the Jeans length at their measured density. This result suggests that most 
of the detected continuum sources are gravitationally unstable and could be collapsing if we 
assume effect only from thermal support. No significant differences are detected between 
prestellar cores and protostellar cores. 

4-1.2. Evolutionary Stages of the Continuum Sources 

The Jeans criterium indicates that most of our SMA sources are gravitationally unstable. 
Here, we examine other observational evidence of their protostellar nature. In Figure 2, we 
compare the positions of infrared sources, i.e., protostars (denoted by orange open circles), 
and T Tauri stars (denoted by white open circles) based on Peterson & Megeath (2008), and 
our 850 fim continuum sources; these comparisons are of similar spatial resolutions (2-5" 
corresponding to 830-2100 AU at Orion). Multi- wavelength data and outflow survey results 
are summarized in column 2-7 of Table 5. Spitzer 8 /im sources were detected toward eight 
out of twelve SMA sources (~67%), six out of twelve SMA sources (~50%) have CO bipolar 
outflow, while three out of twelve SMA sources (~25%) are associated with radio jets (Aso 
et al. 1998; Reipurth et al. 1999; Stank et al. 2002; Williams et al. 2003; Takahashi et al. 
2006; 2008; 2009; Takahashi & Ho 2012). These results suggest that approximately half of 
the SMA sources in this region are protostellar. 

Eight out of 16 Spitzer sources (within the SMA observing region) are associated with 
the SMA continuum sources, which are identified as protostars (Peterson & Megeath 2008; 
Megeath et al. 2012). Six out of 16 Spitzer sources, which are categorized as T- Tauri 
stars have no counterparts in the SMA continuum emission. Two out of 16 Spitzer sources 
identified as protostars (Peterson & Megeath 2008; Megeath et al. 2012) are not associated 
with SMA continuum sources. We did not detect submillimeter continuum emission from 
the T Tauri source as well as two protostars likely due to the mass detection limit of our 
SMA observations. Evolutionary status of the individual sources are summarized in the last 
column of Table 5. Here, the classification between prestellar and protostellar cores is based 
on the detection of either an infrared source, and/or molecular outflow, and/or radio jet. 
The evolutionary status of a protostellar core was determined based on its SED. Since the 
spatial resolution of submillimeter observations are limited, we could not use the original 
definition for identifying Class sources (i.e., L^i/ L sn y )mm < 200 by Andre et al. 1993). 
Instead of this, the classification of Class I (Class 0) sources was determined based on the 
detection (non-detection) of 2.2 jum emission. 
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4.2. Filament Fragmentation 

4-2.1. Hierarchical Structure 

Our SMA continuum observations with 4.5" angular resolution have revealed a chain 
of density enhancements within the OMC-3 filament. The large-scale filamentary structure 
and their fragmentation properties can be a key to understanding the initial conditions of 
embedded clusters such as core spatial distribution, mass distribution, number distribution, 
and evolutionary sequence, all of which may govern the initial core mass function of young 
stellar clusters and initial stellar multiplicity. 

Certain fragmentation length scales within the OMC filaments have been suggested by 
several previous studies (Dutrey et al. 1991; Hanawa et al. 1993; Wiseman & Ho 1998; 
Johnstone & Bally 1999). The largest spacing is shown in the 12 CO image by Maddalena et 
al. (1986), which corresponds to the Orion A and Orion B GMCs, and has a separation of 
~4.7° (35 pc). Within the Orion A cloud, periodical structures have been reported by Dutrey 
et al. (1991) in the C 18 (1-0) image with peak separations about 9' (1.1 pc) toward the 
northern part of the Orion A cloud, while similar structures, with ~10' (1.2 pc) separation, 
have also been reported in the 13 CO (1-0) image at the southern part of the Orion A cloud 
by Hanawa et al. (1993). The emission peaks correspond roughly to significant star-forming 
regions such as OMC-1S, OMC-1, OMC-2, and OMC-3 within the Orion A cloud (John- 
stone & Bally 1999). Single-dish observations in the NH 3 (1, 1) and (2, 2) and single-dish 
(sub) millimeter continuum have shown the 0.1 pc scale dense clumps embedded within the 
OMC filaments (Cesaroni & Wilson 1994; Chini et al. 1997; Lis et al. 1998; Johnstone & 
Bally 1999). Johnstone & Bally (1999) pointed out that these clumps are equally separated 
along the filaments with separation of ~150" (0.31 pc). Moreover, smaller size scale density 
enhancements have been found in the high-angular resolution VLA NH 3 observations of the 
OMC-1 region (Wiseman & Ho 1998). Their mean separation is ?»50" (0.1 pc). 

Our SMA continuum observations revealed fragmentation within the 0.1 pc clumps. 
Nearest neighbor distances for each source are calculated, and listed in the sixth column 
of Table 4. The median separation between nearest neighbors is measured to be ~17" or 
0.035 pc H|. Furthermore, we have performed 1.3 mm continuum observations with the SMA 
toward the OMC -1 north region, and found similar source separation within their filaments 
(«i0.06 pc). The details of the OMC-1 north observations will be presented and discussed in 



3 For the nearest neighbor calculation, we exclude SMM 7 since the separation between SMM 7 and SMM 
8, ~1.4", are much smaller than the other sources. We consider that the separation of the sources are likely 
related to form a binary system, which is the next fragmentation size scale. 
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Teixeira et al. (2013). 

In summary, our SMA observations have resolved the continuum sources with quasi- 
periodical separation as small as ~0.035 pc. A comparisons of the SMA continuum observa- 
tions with previous single-dish observations suggests multi-spatial scale fragmentation within 
the Orion A GMC with sizes scale between 37 pc and 0.035 pc. 



4-2.2. Jeans Analysis 

The fragmentation of molecular clouds is related to gravitational instabilities. Jeans 
analysis can provide a useful way to discuss which scale can grow and produce high-density 
regions. As described in Section 4.1.1, the critical Jeans length is described as equation (4) 
with the assumption of an initially uniform background density (e.g., Jeans 1902; Binney 
& Tremaine 1987). For other geometric configurations, such as cylindrical and sheetlike 
structures, sufficiently large perturbations can also result in gravitational instability (e.g., 
Nakamura, Hanawa & Nakano 1993; Burkert & Hartman 2004). In the case of an infinitely 
long static and cylindrical isothermal cloud, the perturbation wavelength for maximum insta- 
bility (hereafter we simply call this the Jeans length) is described as follows (e.g., Nakamura, 
Hanawa & Nakano 1993; Wiseman & Ho 1998): 



20c s 

In order to understand the physical mechanism of fragmentation processes over a wide 
range of size scale, the measured source separation is plotted as a function of the number 
density of their respective parental molecular clouds/clumps/cores in Figure 5. Note that 
the inclination of the filaments should be taken into account in this discussion. According 
to Genzel & Stutzki (1989), a randomly oriented filament will have a line of sight mean 
(median) inclination angle of 45° (60°). To display this projection effect, (A Q b s = Af rag -sim), 
core separations in the case of % = 45° are plotted in Figure 5 as open symbols. Table 6 
summaries the number densities, structure separations, and gas temperatures of different 
structures. These structures are listed by size: GMCs (>lpc), large-scale clumps (0.5-1.0 
pc), small-scale clumps (0.1-0.5 pc), and dense core (<0.1 pc). Here, we assume that dense 
cores are formed via the fragmentation of parental small-scale clumps, which in turn are 
formed via the fragmentation of parental large-scale clumps, and that these latter structures 
are formed via the fragmentation of parental GMCs. The measured separations were plotted 
as a function of the number density of the parental large scale structure. To further discuss 
the smaller scale fragmentation properties, the separation of the visible binary stars studied 
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in the Orion nebula cluster (ONC; Reipurth et al. 2007) are also plotted assuming that the 
number density of their parental core is similar to that of our SMA continuum source. 

We overplot on Figure 5 the perturbation wavelength, Af rag ., given in equation (2) and 
(4) for various gas temperatures. The observed fragmentation spacings (denoted in Figure 
5 as filled symbols) are roughly following a single power law over almost four orders of 
magnitude, suggesting that the hierarchical fragmentation could be governed by a single 
physical process. Figure 5 indicates that the data is described better by a shallower index, 
-1.4 (dashed line), than that predicted maximum instability size as a function of density of 
an isothermal clouds, -2.0 (solid line). This shallower index can be explained by differences of 
cloud temperatures over different size scales. Parsec-scale GMCs are exposed to ionized gas 
from HII regions, cosmic rays, large-scale turbulence, etc., which work as heating processes. 
Consequently, the mean temperature of parsec scale clouds tends to be relatively high (e.g., 
50-100 K in the ISMs; e.g., Ferriere 2001). Intermediate regions (both large- and small- 
scale clumps) of the OMC molecular cloud (0.1-1.0 pc) is also most likely exposed to the 
strong interstellar radiation field from nearby star forming clusters, e.g., strong heating from 
the exterior, corresponding to 1000 times stronger interstellar radiation field as suggested 
toward a few small-scale clumps in the OMC-2/3 region from the CO isotopes observations 
with Monte Carlo line radiative transfer models (e.g., J0rgensen et al. 2006). This strong 
radiation is able to produce the higher gas temperature (>30 K) at the outer most part of the 
envelope gas. Finally, deeply embedded dense cores within filaments are shielded from these 
heating processes. Thus, the innermost parts of clouds have low temperatures (e.g., ^20 K; 
Cesaroni & Wilson 1994; Wiseman & Ho 1998; Tatematsu et al. 2008). This decrease in 
temperature with size scale, means the fragmentation length will also decrease more sharply 
than the isothermal case. This leads to a more shallow power law in Figure 5. 

Moreover, the dissipation of the magnetic flux on the smaller scale can also be responsi- 
ble for the observed shallower power law index. Since large-scale diffuse structures are always 
exposed to cosmic ray and/or ionized gas from the surrounding HII region, the fractional 
ionization in large-scale molecular clouds is higher than in the small scale dense gas core. 
Therefore, ambipolar diffusion will be more efficient in the dense cores, leading to less mag- 
netic support, and a smaller fragmentation length (e.g., Nakano & Tademaru 1972). Finally, 
turbulence, which is often observed in GMC size scales (Larson 1981; Solomon et al. 1987), 
can also be responsible for the observed shallower power law index. Detailed discussion for 
the OMC case will be described in the Section 4.2.3. 

Observational results suggest that the fragmentation spacings, especially at the scales 
of clumps (0.1-1 pc), are roughly consistent with the maximum instability size derived from 
either the initially uniform density cloud (orange solid line in Figure 5) or cylindrical clouds 
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(blue solid line in Figure 5) assuming increasing temperatures with size scale. This indicates 
that the fragmentation may be occurring through the thermal fragmentation process. The 
SMA detected continuum sources (denoted in Figure 5 by triangle shaped data points) are 
embedded within a filament. Hence the core separation should be around the maximum 
instability size of the cylindrical clouds model (denoted in Figure 5 by blue line). However, 
the observed core separation is a factor of 2.0-7.3 smaller than the thermal fragmentation 
spacing with T gas =20 K, n Ha = 8.8 xlO 4 - l.lxlO 6 cm~ 3 (Table 6), and i=90°. In the next 
subsection, we explore different physical processes that can affect the fragmentation size 
scale. 



4-2.3. What Mechanisms affect the Observed Core Separation? 

From our SMA observational result, the observed fragmentation (Af rag =0.035 pc) is 
smaller than the Jeans length of the clump. Projection effects are not sufficient to explain 
this difference: Af rag = A obs /sim < Aj for 45°>i>90°. In order to further discuss the frag- 
mentation properties within the Orion A GMC, we summarize different possible mechanisms 
that can affect the fragmentation size scale. 

Turbulence: In this case, the cloud fragmentation length is described by equations 
(2) and (4) with the effective sound speed c cS instead of the sound speed c s . The c eS 
is defined as c^ ff = A ${£f > wnere At> bs is the observed velocity dispersion and is given 
by At>Q bs = Avf + Av„ t , where Av t and Av nt correspond to the thermal and non-thermal 
velocity dispersion, respectively. Turbulence will therefore produce a fragmentation scale 
larger than the Jeans length. Figure 5 and Table 6 present that observed clump separations 
are roughly consistent, or smaller, than the thermal fragmentation lengths for scales of a few 
pc. These results suggest that turbulence has already decayed for size scales <1 pc. For scales 
of a few pc or larger, turbulent pressure may in fact be dominant. For example, Jackson et 
al. (2010) studied a filament that spans 80 pc and conclude that its observed fragmentation 
scale (~5 pc) can only be described by an isothermal filamentary cloud determined by 
turbulent pressure, or equivalently, by a self-gravitating incompressible fluid cylinder. Our 
data indicates that the OMC 3 filament fragmentation cannot be explained by this scenario. 

Magnetic field and cloud rotation: Other physical processes that might affect the 
fragmentation length of filaments are helical magnetic field and cloud rotation (Fiege & 
Pudritz 2000; Matsumoto, Nakamura, & Hanawa 1994). In fact, Matthews et al. (2001) 
use helical field geometry of magnetic field to interpret their JCMT/SCUBA polarization 
observations toward OMC-3. Their data indicate that the magnetic field in the northern 
part of OMC-3 is toroidal. Hanawa et al. (1993) has performed magnetohydorodynamical 
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simulations for filament fragmentation, consisting of an isothermal rotating cloud + helical 
magnetic field, and found a dispersion relation (eq 11, Hanawa et al. 1993). We can write 
the fragmentation wavelength for the toroidal magnetic field case as: 

frag 0.72[(1 + a/2 + 0)V3 - 0.6] 1 ' 

where H represents a length scale of the parental clump, and given by 

where p c is the peak density of the parental clump, and the effect from magnetic field and 
rotation are described by a and f3 as: 

B 2 a 2tt 2 H 2 

a = ^r— P = — — ( 7 ) 



8n p c c 2 s ' 



Q is the rate of rotation and B is the magnetic field strength. Equation (5) shows that if 
(a + (3)>l, the maximum instability size will become smaller than the thermal fragmentation 
length. If a + ft « 1, the effects from the magnetic field and rotation are negligible. 

Global collapse of the parental cloud Recent theoretical studies present how cloud 
structures and global collapse significantly affect local fragmentation properties (e.g., Burket 
& Hartmann 2004; Pon et al. 2011; Pon et al. 2012). If the original parental cloud is not 
infinitely long or spherical, gravitational focusing results in enhanced concentrations of mass 
inducing global collapse unless magnetic fields, cloud rotation, and/or external pressure are 
dominated (Burket & Hartmann 2004). As we can see from the JCMT/SCUBA 850 /im 
continuum observations (figure 2, right panel), the parental filamentary clump (e.g., 1.5 
Jy beam -1 denoted by thick contour) shows a relatively well defined edge, so according to 
Burket & Hartmann's model, global collapse could potentially be underway in this clump. A 
recent paper by Pon et al (2011) pointed out that filamentary geometry is the most favorable 
situation for the smallest perturbation to grow before global collapse overwhelms the cloud. 
Their following paper (Pon et al. 2012) suggested that filaments with large aspect ratio 
would take more time to collapse. For example, the global collapse time-scale for a cloud 
aspect ratio of 10 (similar case of the observed clumps) is roughly seven times slower than the 
spherical free-fall time. This can explain how local fragmentation at the size scale of dense 
cores and global collapse of their parental small-scale clump could be concomitant. Global 
collapse will eventually shorten the spacing between the fragmented cores, and this scenario 
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could therefore potentially explain the observed shorter fragmentation length compared to 
Jeans length. 

In summary, there are three physical mechanisms that can increase (turbulence) or de- 
crease (magnetic fields and rotation, or global collapse) the fragmentation length relative to 
the Jeans length. The fragmentation lengths obtained from the data presented in this paper 
are consistent (0.3 pc -lpc scale) or smaller (< O.lpc scale) than the Jeans fragmentation 
length, suggesting that turbulence does not play a significant role in the fragmentation pro- 
cess for these size scales. Note that some of their core are larger than their associated Jeans 
length, therefor, they could fragment further. We discussed how the magnetic field and cloud 
rotation, or global collapse may decrease the observed fragmentation length. However, at the 
moment we are unable to quantitatively constrain a + (3, or inflow from the global collapse 
for our case - only future direct comparisons between higher- dynamic range kinematic data 
and more realistic models can properly address these issues. 

Core Fragmentation: In addition to the clump large-scale fragmentation properties 
discussed above, we would like to explore further the smaller core-scale fragmentation prop- 
erties, which is related to binary star formation (<1000 AU). These binary stars observed 
in the Orion Nebula Cluster (ONC) have a wide range of separations from 0.0027-0.027 pc 
(67.5-675 AU; Reipurth et al. 2007). The plotted results in Figure 5 suggest that wide binary 
systems (>200 AU) might follow the same power law index extrapolated from the observed 
large size scale cloud fragmentation, i.e., these wide binaries could be formed through the 
next level of hierarchical process. This lends support to the theoretically predicted binary 
formation scenario through the hierarchical fragmentation of the parental core occurs in the 
isothermal collapse phase (e.g., Goodwin et al. 2007; Machida et al. 2008). Close binary 
systems (<100 AU) might be formed through other mechanisms; hydrodynamic simulations 
suggest that angular momentum within the rotating and collapsing hydrostatic first cores is 
a crucial parameter to determine the later fragmentation properties (e.g., Boss 1989; Bon- 
nell and Bate 1994; Bate 1998; Saigo & Tomisake 2006; Machida et al. 2008). Furthermore, 
star-disk encounters could also play an important role in re-distributing angular momentum 
in protoplanetary disks and forming binary systems (e.g., Larson 2002; Pfalzner 2004; 2005). 



4-2.4- Fragmentation Timescale 

One interesting issue related to filament fragmentation is whether the fragmentation 
occurs simultaneously in the entire filament, or instead propagates along it producing se- 
quential star formation along the filament. With the SMA continuum observations, we found 
the following interesting tendencies among the sources; (i) a chain of continuum sources is 
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detected in the northern part of the OMC-3 filament from SMM 2 to SMM 8 with a similar 
spatial separation, (ii) a relatively extended fluffy source (elongated structure), SMM 5, is 
located at the center of the chain of the continuum sources, and this fluffy structure seems 
to be still in the process of core formation (or fragmentation). No star formation signatures 
such as outflow/jet association or infrared source are found, which supports the idea that 
SMM 5 is probably the youngest star formation site within the filament. Why does SMM 5, 
which is located in the middle of the filament, have the fluffy structure? A fragmentation 
process which propagates along the filament from each end inwards can be one possibility 
to explain this. 

Previous numerical calculations give us a hint that the finite geometry of gas clouds 
promotes concentrations of material at the ends of the filament, at the beginning of the 
filament fragmentation. Hanawa, Yamamoto, & Hirahara (1994) performed one-dimensional 
hydrodynamic simulations of a finite long filamentary cloud. These simulations suggest that 
dense cores are produced with a semi regular interval in space and time from one edge to the 
other. The gas near the edge is attracted inwards by gravity and the accumulation of the 
gas makes a dense core near the edge. When the dense core grows in mass up to a certain 
amount, it gathers gas from the other direction. Accordingly the dense core becomes isolated 
from the main cloud and the parent filamentary cloud has a new edge. Via the same process, 
the next accumulation of gas at the new edge makes a new dense core. The fragmentation 
is considered to propagate with approximately the effective sound speed (1.23c e fr derived 
from the simulation by Hanawa Yamamoto & Hirahara 1994; l.lc e ff derived from a linear 
stability analysis by Nakamura, Hanawa, & Nakano 1993). Signatures of the fragmentation 
propagating down the filament is actually observed in the Taurus Molecular Cloud -1 region 
(Hanawa, Yamamoto & Hirahara 1994). 

Our SMA results also show that the chain of continuum sources, from SMM 2 to 
SMM 8, are embedded within the most northern part of the OMC-3 filaments where both 
ends show edge-like density contrast (finite filamentary clouds with two edges; Figure 2). 
The mean value of the observed velocity width in FWHM is approximately 0.88 km s _1 
in this region (A mean velocity width derived from NH 3 observations toward the norther 
OMC cloud by Cesaroni & Wilson 1994) and the effective sound speed is estimated to be 

c c fj~-^/^^~0.37 km s -1 . The mean projected distance from the sources at the two edges 
of the filament (i.e., SMM 2 and SMM 8) to SMM 5 is measured to be 0.15 pc. There- 
fore, the time scale of the fragmentation propagation is estimated to be approximately 
0.15 pc/0.37 km s^S.OxlO 5 yr. The detected continuum source in the filament have sim- 
ilar to the protostellar phase (Peterson & Megeath 2008; Megeath et al. 2012). The typical 
evolutionary time scale of such sources are 10 4 -10 5 yr (Beichman et al. 1986; Shu et al. 
1977). Estimated propagation time scale of the filament fragmentation is similar to, or 
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slightly longer than, the star formation time-scales of the detected sources. This may sug- 
gest that the protostars formed within the fragmented cores, where star formation began 
almost at the same time as the clump fragmentation. More complete fragmentation time 
scale analysis (including other parts of the OMC filaments) will be presented in following 
papers. 

5. CONCLUSION & FUTURE PROSPECTS 

High-resolution 850 /zm continuum observations have been carried out using the SMA 
toward the OMC-3 northern filament with a mosaic of 85 fields covering approximately 
6'.5x2'.0 (0.80x0.25 pc). The data were supplemented by multi-wavelength data sets for 
investigating evolutionary stage of each source. The main results are summarized as follows: 

1. With the SMA 850 /im observations, 12 dust continuum sources have been identified 
and spatially resolved with an H 2 mass between 0.3 -5.7 M and the source size be- 
tween 1400- 8200 AU. We found 850 jum continuum counterpart to the previously 
known ten sources and present the discovery of two new additional sources. The 
detected sources have a variety of structures and evolutionary stages, ranging from 
prestellar phase to the Class 0/1 phases. The Jeans analysis suggests that all the 
detected sources are most likely gravitationally unstable. Multi-wavelength data sets 
confirmed that 6/12 (50 %) continuum sources have molecular outflows, 8/12 (67 %) 
continuum sources are associated with infrared sources, and 3/12 (25 %) continuum 
sources are associated with ionized jets. These results suggest that the OMC-3 is an 
active protostar forming region. 

2. Hierarchical structure with quasi-periodical separations on different size-scales is sug- 
gested with a size of ~4.7°/35 pc down to 17"/0.035 pc, implying a multi-size scale 
fragmentation process. Fragmentation spacings in clumps (0.1-1 pc) are roughly con- 
sistent with the thermal fragmentation lengths in each structure. However, shorter 
spacings than those expected from thermal fragmentation have been detected in the 
SMA data (^0.04 pc). This could be explained by either an helical magnetic field, 
parental cloud rotation, or/and large-scale filament collapse. 

3. Possible evidence for sequential fragmentation as local conditions evolve is suggested in 
the northern part of the OMC-3 filament. The fragmentation propagation timescale is 
estimated to be r fl . ag .~3.0xl0 5 yr. This time scale is roughly consistent with (or slightly 
longer than) the life time of the detected protostellar sources, which is 10 4 -10 5 yr. 
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These results suggest that the protostars formed within the fragmented cores, and 
that star formation began at almost the same time as the core fragmentation. 

Our interferometric observations successfully resolved the detailed gas structure asso- 
ciated with the individual star formation sites within the filamentary structure. However, 
present observations are still limited by mass sensitivity (>0.15 M ) and suffer from the 
large missing flux of the filament (~16% recovering flux). Due to this, we are biased to de- 
tect relatively bright and compact structures (i.e., mostly massive envelope associated with 
the pre-stellar and protostellar cores). In order to detect less bright sources (e.g., emission 
from Class II sources, brown-dwarf mass sources) and to obtain a wider evolutionary range 
from the sample, even higher sensitivity observations are required. In addition, a larger 
sample is crucial to investigate the evolutionary tendencies along the entire OMC filaments. 
Furthermore, polarization observations with high sensitivity and high-spatial resolutions are 
crucial for elucidating the magnetic field structures, for measuring the field strength at the 
early evolutionary stages, and for verifying the role of the magnetic field in the filament 
fragmentation processes. These issues are future prospects with extended data sets using 
the SMA, EVLA, and ALMA. 
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Fig. 1. — Left panel: An 850 /im dust continuum image taken with the JCMT/SCUBA, 
referred from Johnstone & Bally (1999). Right panel: SMA 85 points mosaic. The back- 
ground image is same as in the left panel. Open circles show the SMA primary beam at 
351 GHz. Crosses in the image show the positions of previously detected 1.3 mm continuum 
sources with a single-dish telescope (Chini et al. 1997). 
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Fig. 2. — Left panel: An 850 /im continuum mosaic image taken with the SMA. White 
contours show the 5a signal levels (i.e., 100 mJy beam -1 ). Orange open circles and 
white open circles present the positions of protostars and T Tauri stars, respectively 
(Peterson & Megeath 2008). Right panel: The SMA 850 /im continuum mosaic image (color) 
overlaid with the 850 fim continuum image taken with the JCMT/SCUBA (contours referred 
from Johnstone & Bally 1999). The contour levels start at 0.3 Jy beam -1 with a interval of 
0.6 Jy beam -1 . The orange and white circles are the same as in the left panel. Black filled 
circles in the bottom right corner show the JCMT and SMA beam sizes of 14" and 4". 5, re- 
spectively. Thick contour (contour level of 1.5 Jy beam -1 ) presents a relatively well defined 
edge of the small-scale clump discussed in Section 4.2.3. 
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Fig. 3. — Zoomed-in images of SMA detected continuum sources. The contour levels of 
these images (except the SMM 7-9 panel) start at 5cr level (i.e., 100 mJy beam -1 ) with a 
interval of 5cr. Contour levels of SMM 7-9 show 5, 10, 15, 30, 45, 60, ... 210 a. In these data 
SMM 7 is not resolved, Takahashi et al. (2009) shows the source resolved and accurately 
positioned. Positions of 8/xm Spitzer sources identified by Peterson & Megeath (2008) are 
denoted by orange crosses. The bars in the bottom right corner in each panel show the 
scale of 2100 AU, which is the SMA beam size in FWHM. Color wedge denoted at the SMM 
4 panel shows the signal levels of all the individual panels. The filled circles in bottom 
right corner of the central panel presents the SMA and JCMT beam sizes of 4". 5 and 14", 
respectively. 
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Fig. 4. — (A) Histogram of source mass. The SMA detection limit of 1.5 M was estimated 
from equation 1 with the following assumptions: total flux of 100 mJy (5 a), Td ust =20 K,and 
(3=1.6. (B) Histogram of equivalent circular source size, =\J D^ i&) + -D^ in , where -D ma j. 
(-Dmin.) is size of the major (minor) axis of the source in AU. (C) Filled and open circles 
show the source sizes as a function of number density for protostellar and prestellar sources. 
Dotted, solid, and dashed lines show the critical Jeans length as a function of number density 
with the gas temperature of T gas =10, 20, and 30 K, respectively. 



-30 - 




0.0001 0.001 0.01 0.1 1 10 100 

Separation (PC) 



Fig. 5. — Observed separations of clouds/large-scale clumps/small-scale clumps/cores in 
parsec as a function of the mean hydrogen number density of the parental cloud in cm -3 . 
The filled symbols show the separations with an assumption of % = 90° (filament in the plane 
of the sky), and the open symbols show the separation with an assumption of % = 45°. The 
plotted values are summarized in second and fourth columns of Table 6 with reference papers. 
Data from GMC (Sakamoto et al. 1994; Wilson et al. 2005), large-scale clumps (Dutrey et 
al. 1991; 1993; Hanawa et al. 1993; Johnstone & Bally 1999), small-scale clumps (Cesaroni 
& Wilson 1994; Johnstone & Bally 1999), dense cores (This work), and binaries (Reipurth et 
al. 2007) are denoted by diamonds, circles, squares, upright triangles, and inverted triangles, 
respectively. The dashed line shows a reference power law index of (nH 2 ocA -L4 ). The color 
lines show the expected maximum instability size as a function of number density derived 
from equations (4) and (5): uniform background density with a gas temperature of 20 K 
(orange), cylindrical cloud with a gas temperature of 20 K (blue line), 50K (green), and 100 
K (pink). 
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Table I. SMA Observational Logs 



Date 


Configuration Number of Antennas 


i/ c (USB)/i/ c (LSB) 
(GHz) 


T225 


Tsys[DSB] 

(K) 


oJ0423-013 
°351GHz 

(Jy) 


oJ0530+135 
°351GHz 

(Jy) 


2008 Jun 2 


compact 7 


356.47/346.47 


0.04-0.1 


120-460 


1.61 


1.58 


2008 Aug 29 


subcompact 8 


356.55/346.55 


0.04 


100-200 


1.05 


1.43 


2008 Oct 21 


compact 8 


356.47/346.47 


0.08-0.1 


240-450 


1.06 


1.82 



Table 2. SMA Continuum Observational Parameters 



Parameter Value 

Reference position (J2000.0) a=5 h 35 m 23.48 s , <5=-5°01'32.20" 

Configurations 2 X compact and lx subcompact 

Primary beam HPBW (arcsec) 36 

Grid spacing (arcsec) 18 (P.A. 45°) 

Synthesized Beam HPBW (arcsec) 4.5x4.5 

Equivalent frequency (GHz) 351 

Total continuum bandwidth USB+LSB (GHz) 3.2 

Projected base line range (kA) 9.7—88 

Maximum detectable structure (arcsec) a 17 

Gain calibrators 0423-013 and 0530+135 

Bandpass calibrator 3C 273 and 3C 279 

Flux calibrators Titan and Neptune 

RMS noise level (mjy beam -1 ) ~20 



a Our observations were insensitive to more extended emission than this size-scale 
structure at the 10% level (Wilner & Welch 1994). 



Table 3. SMA 850 /im Continuum Sources 



ID 


R.A. a 


Dccl. a 


Total Flux 


Peak Flux f 


Deconvolved Size 


P.A. 


C(SMA/JCMT) b 


Other mm/submm source names c 




(J2000.0) 


(J2000.0) 


(mjy) 


(mjy beam -1 ) 


(arcsec) 


(degree) 




OMC3-SMM 1 


5 35 16.03 


-5 00 03.07 


354±63 


148±19 


8.5±1.0x2.3±l.l 


125±5 


0.23 


CSO 3 


OMC3-SMM 2 


5 35 18.03 


-5 00 17.77 


1070±42 


811±20 


2.8±0.2x2.2±0.2 


107±13 


0.33 


MMS 1, CSO 5 


OMC3-SMM 3 


5 35 18.30 


-5 00 33.01 


1001±43 


736±20 


3.9±0.2xl.0±0.4 


81±3 


0.36 


MMS 2, CSO 6 


OMC3-SMM 4 


5 35 18.98 


-5 00 51.63 


1009±52 


220±21 


9.5±1.9x7.6±1.7 


-36±38 


0.31 


MMS 3, CSO 7 


OMC3-SMM 5 


5 35 20.88 


-5 00 56.25 


2572±65 


393±30 


19±1.6x5.2±0.6 


-55±3 


0.49 


MMS 4, CSO 8 


OMC3-SMM 6 


5 35 22.47 


-5 01 14.43 


1261±47 


810±20 


3.9±0.2x2.8±0.2 


100±6 


0.37 


MMS 5, CSO 9 


OMC3-SMM 7 


5 35 23.50 


-5 01 29.58 


1600 d 


320 d 


5.5x3.8 d 


-15 


_C 


MMS 6-NE 


OMC3-SMM 8 


5 35 23.42 


-5 01 30.35 


6152±44 


4400±20 


3.9±0.03xl.4±0.1 


73±1 


0.86 


MMS 6, CSO 10 


OMC3-SMM 9 


5 35 23.49 


-5 01 39.21 


874±60 


3900±19 


5.2±0.3x4.8±0.3 


150±28 






OMC3-SMM 10 


5 35 28.20 


-5 03 40.65 


220±53 


118±20 


5.1±l.lx3.2±1.2 


140±23 


0.28 


MMS 7-NE 


OMC3-SMM 11 


5 35 26.56 


-5 03 55.04 


1162±41 


921±20 


3.2±0.2xl.0±0.4 


52±4 


0.56 


MMS 7, CSO 12 


OMC3-SMM 12 


5 35 25.44 


-5 03 53.22 


290±64 


117±19 


7.3±1.2x3.8±1.2 


80±12 


0.32 





oo 

a The source positions were determined from the non-self calibration images with emission peak positions. The positional difference between the non self-calibrated ^ 
images and self-calibrated images is «0.1". 

b The flux concentrated ratio is calculated from the total flux measured from the SMA image divided by the JCMT flux within 14" beam. SMA image are convolved 
with the JCMT beam size of 14". The uncertainties of the derived ratios are ±0.04 (in absolute value). 

c REFERENCE. — MMS sources by Chini ct al. (1997); CSO sources by Lis ct al. (1998); MMS 7-NE by Takahashi ct al. (2006); MMS 6-NE by Takahashi et al. (2009) 

d Flux information and deconvolution size of SMM 6-NE are referred from Takahashi ct al. (2009) 

e Single-dish data taken with the JCMT were not able to fit with the Gaussian components due to the separation between sources compared with the large spatial 
resolution. 
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Table 4. Physical Parameters 



ID 


M H2 a 


Nh 2 * 




Source Size b 


NN distance 




(Af©) 


(10 23 cm" 2 ) 


(10 7 cm- 3 ) 


(AU) 


(arcsec) 


SMM 1 


0.5±0.1 


0.9±0.4 


0.3±0.2 


3520x950 


33.3 


SMM 2 


1.6±0.3 


17.0±4.4 


15.0±4.5 


1160x910 


15.8 


SMM 3 


1.5±0.3 


13.0±3.1 


9.6±2.8 


1620x410 


15.8 


SMM 4 


1.5±0.3 


1.4±0.7 


0.4±0.2 


3930x3150 


21.2 


SMM 5 


3.9±0.8 


1.4±0.4 


0.2±0.1 


7870x2150 


28.7 


SMM 6 


1.9±0.4 


11.0±2.5 


7.1±1.8 


1620x1160 


21.3 


SMM 7 


2.4±0.5 


7.3±1.5 


3.3±0.7 


2280x1570 




SMM 8 


5.7±l.l d 


44.0±8.9 d 


33±6.6 d 


1620x580 


8.9 


SMM 9 


1.3±0.3 


3.6±0.9 


1.5±0.4 


2150x1990 


8.9 


SMM 10 


0.3±0.1 


1.2±0.8 


0.6±0.5 


2110x1330 


28.4 


SMM 11 


1.8±0.4 


21.0±5.3 


19.0±5.8 


1330x410 


16.8 


SMM 12 


0.4±0.1 


0.9±0.4 


0.3±0.2 


3020x1570 


16.8 



a Measurement errors of deconvolved size and total flux described in Table 3 arc taken 
into account for error estimations. 

b Tlie error for source sizes and position angles are summarized in Table 3 column 
5-7, respectively. 

c The nearest neighbor distance for each source (Section 4.2.1). SMM 7 was excluded 
from the nearest neighbor calculation since the separation between SMM 7 and SMM 
8 is much smaller than angular resolution of this paper (4". 5). We consider that the 
separation between SMM 7 and SMM 8 are related to formation of a binary system. 

d /3=0.93 derived by Takahashi et al. (2009) was adopted to estimate the parameters. 
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Table 5. Comparison of Multi-wavelength Data Sets toward Detected Continuum Sources 



ID 


CO outflow H2 jet 


NIR nebula 


3.6 cm emission X-ray 


Infrared Source Detection 


Evolutionary Status a 


SMM 1 








2.2, 3.8, 4.5, 5.8, 8, 24 9 ' J0 


Protostellar core 


SMM 2 










Prestellar core 


SMM 3 


Yes 1 ' 2 Flow B 5 


Yes 2 


VLA I 6 - 7 Yes s 


2.2, 3.8, 4.5, 5.8, 8.0, 24 9 ' 10 


Protostellar core 


SMM 4 






Yes« 


4.5, 5.8, 8, 24 9 - 10 


Protostellar core 


SMM 5 










Prestellar core 


SMM 6 


Yes 1 ' 2 Flow C 5 


Yes 2 




3.6, 4.5, 5.8, 8, 24 9 ' 10 


Protostellar core 


SMM 7 




Yes 2 


VLA 3 6 


2.2, 3.8, 4.5, 5.8, 8.0, 24 9 ' 10 


Protostellar core 


SMM 8 


Yes s 






2 4 9,io ? b 


Protostellar core 


SMM 9 




Yes 2 




2.2, 3.6, 4.5, 5.8, 8, 24 9 < 10 


Protostellar core 


SMM 10 


Yes^ 


Yes 2 




3.6, 4.5, 5.8, 8.0, 24 9 ' 10 


Protostellar core 


SMM 11 


Yes 1 ' 2 '* Flow F' 5 


Yes 2 


VLA 4 6 


2.2, 3.6, 4.5, 5.8, 8, 24 9 ' 10 


Protostellar core 


SMM 12 










Prestellar core 



REFERENCE. — 1: Willams et al. (2003), 2: Takahashi et al. (2008), 3:Takahashi & Ho (2012), 4: Takahashi et al. (2006), 5: Yu et al. 
(1997), 6: Reipurth et al. (1999), 7: Tsujimoto et al. (2004), 8: Tsuboi et al. (2000), 9: Peterson & Megeath et al. (2008), 10: Mcgcath et 
al. (2012) 

a Evolutionary status of each source was referred from Megeath et al. (2012). 
b Specially unresolved source at 24 /im. 



Table 6. Hierarchical Structure of the OMC filaments 





riH 2 


Tgas 


Measured Separation 


Corresponding Size Scale 




(cm^ 3 ) 


(K) 


(arcminute) 




ISMs 


50-100 J 


50-100 J 




ISMs observed in HI 21 cm observations 


GMCs 


200-3000 2 




282 s 


Ori A and Ori B (Separation between Orion BN/KL and NGC 2024) 


Large-scale Clumps 


9.0e03-3.0e04^ 




9-10 5 ' 6 ' 7 


Significant star forming regions such as OMCl,2,3, etc. 


Small-scale Clumps 


8.8e04-l.le06 5 '* 


15-28 8 


2.5 s 


identified with single-dish NH3 and continuum observations 


Dense Cores 


3.2c06-3.3c08 9 




0.15-0.55 9 


Individual star forming sites detected with the SMA observations 


Binaries 






0.0025-0.025 w 


Binary separation range measured in ONC (Ha observations). 



REFERENCE. — 1: Ferriere (2001), 2: Sakamoto ct al. (1994), 3: Wilson ct al. (2005), 4: Dutrey ct al. (1993), 5: Johnstone & Bally (1999), 6: 
Dutrey et al. (1991), 7: Hanawa et al. (1993), 8: Ccsaroni & Wilson (1994), 9: This work, 10: Reipurth ct al. (2007) 



